Metabarcoding and Metagenomics

The marine biodiversity inhabiting rocky shores in the Colombian Pacific remains largely undocumented, primarily due to geographic isolation, logistical challenges, and socio-political constraints. To address the existing knowledge gap, we conducted an expedition to enhance baseline biodiversity knowledge in rocky shores by integrating multiple complementary approaches, including visual censuses, specimen collection with morphological identification, environmental DNA (eDNA) metabarcoding and DNA barcodes. eDNA samples were collected at four coastal sites adjacent to rocky substrates, along with biological specimens obtained from fourteen locations through SCUBA diving at depths ranging from 1 to 25 meters. Tissue samples were subjected to genomic DNA isolation, followed by the generation and validation of cytochrome c oxidase subunit I (COI) barcode sequences, which were subsequently corroborated through taxonomic assessment to ensure accurate species identification. eDNA metabarcoding analyses yielded over 7 million high-quality sequence reads. Although taxonomic resolution at the species level was constrained by the limited completeness of reference sequence databases, a total of 106 species and 83 families were successfully identified, predominantly within the classes Actinopteri, Chondrichthyes, and marine mammals. From the 769 specimens obtained we generated 871 sequences, including 414 validated COI barcodes representing 76 species across 64 families. The integration of DNA barcoding and eDNA approaches resulted in over 1,400 taxonomic detections spanning five phyla, with only six species shared between methodologies. Richness and diversity varied among sites, and revealed significant differences along the coastline between Jurado and Cupica Gulf. All sequences were deposited in BOLDsystems database under the CCBIO project and were visualized through OBIS and GBIF databases. These findings provide the first molecular-based baseline for rocky shore biodiversity in the Colombian Pacific, highlighting the value of integrative approaches for monitoring and conservation.

Molecular identification of cetacean bone remains can be limited by DNA degradation and the presence of PCR inhibitors. Here, we present an optimized DNA extraction protocol based on a total demineralization method for environmentally exposed cetacean bones. The protocol uses 100 mg of bone powder, 24 h digestion with EDTA, N-lauroylsarcosine, and proteinase K, followed by a modified silica-column purification. Nine environmentally degraded bone samples representing eight individuals were processed. DNA concentrations ranged from 7.3 to 57.1 ng/uL (mean SD = 25.91- 13.91 ng/uL). The mitochondrial cytochrome b gene was successfully amplified from all samples using conventional PCR, and five samples (55.6%) yielded sequences suitable for downstream analysis. BLASTn identified Balaenoptera physalus as the closest database match for all recovered sequences, and phylogenetic analysis further supported their association with B. physalus reference sequences. These results demonstrate that the proposed protocol provides a practical approach for recovering amplifiable and molecularly informative mitochondrial DNA from environmentally degraded cetacean bone material, facilitating molecular identification from challenging skeletal remains.

Bottom trawling is the most widespread human disturbance on the seafloor, known to affect seabed structure and benthic biota. Yet, its impact on benthic communities remains documented mainly for taxa identifiable through morphology, leaving microscopic benthic communities largely overlooked. The Gulf of Lion continental shelf has been trawled for decades by French and Spanish fleets, to the point that encountering untrawled reference areas became elusive. Here, we examined prokaryote, protist and Animalia benthic communities at eight stations, using eDNA metabarcoding of core sediment samples, with five primer-pairs. Benthic communities were cross-referenced with contrasted trawling impacts, integrated across three periods: over one month, three years or ten years preceding the sampling. The community composition was primarily structured by ecological habitats, yet, trawling pressure showed association with prokaryotes and Animalia diversity. Diversity was higher at sites trawled more intensively in the preceding month, whereas long-term trawling was associated with a significant erosion of Animalia alpha diversity. Unexpectedly, communities did not display consistent vertical structure and the relative contributions of trawling and hydrodynamic reworking remain unresolved due to the absence of untrawled sites. Overall, our results show that trawling effects extend to prokaryotic compartments rarely considered in impact assessments, providing a proof of concept for eDNA-based assessment, and underline a pervasive long-term effect on Animalia communities. This study also highlights the importance of the timescale used to quantify fishing pressure and underscores the complexity of impact assessment when monitoring begins too late to establish a clear baseline, particularly when no reference site remains.

DNA barcoding involves the recovery of a DNA sequence for a target gene region from its source specimen. This process gains complexity when multiple sequences are recovered from a specimen, as is often the case when data are generated by high-throughput sequencers. This diversity can reflect both methodological artifacts (e.g., chimeras, PCR errors, sequencing errors, tag jumps) and real template diversity in the DNA extract (e.g., contamination, endosymbionts, NUMTs, parasites). To support analysis of the sequence data from three million specimens annually, the Centre for Biodiversity Genomics (CBG) has developed BIP, the Barcode Inference Pipeline. Compatible with all sequencing platforms, BIP processes .fastq files and returns both target DNA barcodes and non-target sequences. To generate results, BIP implements quality and size filtration, demultiplexing, primer trimming, chimera scanning, sequence error correction, OTU delineation, and sequence identification. When analysis targets the cytochrome c oxidase 1 (COI) barcode region, BIP also assigns each OTU to a known BIN or identifies its nearest neighbour BIN. As final output, BIP returns summary files ready for upload to BOLD or for other downstream analyses. They include a taxonomic assignment for each OTU, generated by comparison with a DNA barcode reference library. We describe BIPs flexibility and structure, then demonstrate its functionality by analyzing COI sequence data from 100K specimens. Because of its capacity to disentangle target and non-target sequences, BIP outperforms an alternative software package, ONTbarcoder, in several important ways. To ease access, installation, and functionality, BIP is provided as a Docker container (github.com/cbg-innov/BIP).

Environmental DNA (eDNA) has emerged as a transformative tool for biodiversity monitoring, offering a non-invasive means to detect and identify species across various ecosystems. However, despite its potential, current eDNA methodologies have significant limitations, including limited sensitivity, restriction to a specific taxonomic group, lengthy processing times, the need for specialized laboratory conditions and equipment, as well as computational bottlenecks in equipment, connectivity, and algorithms. Overcoming these limitations is crucial for the real-time monitoring of a high number of taxa in field-based scenarios. Here we present an integrated eDNA workflow from extraction to data analysis that mitigates several of these challenges through several technical innovations. Our approach adapts the Tagsteady protocol to Nanopore sequencing, integrates multiplex PCR for simultaneous plant and animal detection, and introduces RTBioScan, a real-time analysis pipeline employing hierarchical basecalling and OTU-level consensus generation to achieve high-accuracy taxonomic assignments during sequencing. Using a single set of universal primers, we integrate metabarcoding and barcoding into a cohesive protocol, utilizing universal markers for both animals and plants (cytochrome oxidase 1 region, COI, and Internal Transcribed Spacer 2, ITS2). This comprehensive approach can detect and identify a broad spectrum of species, from individual organisms to entire communities in a single workflow.

Background: Amplicon sequencing of the 16S rRNA gene is widely used to characterise microbial communities, but the performance of commonly used clustering and denoising pipelines has not been systematically evaluated for highly diverse environmental datasets. We therefore benchmarked four established clustering and denoising pipelines, VSEARCH cluster_size, UNOISE as implemented in VSEARCH, Swarm, and DADA2, using simulated microbial communities with known ground-truth compositions and evaluated how performance varied with species richness, abundance unevenness, sequencing depth, and minimum abundance threshold. Results: Differences among pipelines were small at low species richness but became greater as species richness increased. High average cluster purity did not necessarily correspond to accurate species-level recovery, as species could be split across multiple clusters or recovered incompletely. Among the pipelines, UNOISE may be preferable when high species representation and cluster purity are prioritised, while also showing strong reconstruction of relative abundance profiles, albeit with extensive species splitting and many unclustered reads. DADA2 showed similarly strong reconstruction of relative abundance profiles and less species splitting but represented fewer species and had lower cluster purity at higher richness. Swarm provided a balanced compromise between species representation and splitting, whereas cluster_size may be useful when limiting species splitting is a priority, despite weaker reconstruction of relative abundance profiles. Increasing the minimum abundance threshold reduced species splitting but also reduced the number of clusters and perfect clusters and, at higher thresholds, decreased concordance with ground-truth relative abundance profiles. These effects were more pronounced at lower sequencing depths. Analysis of seafloor sediment samples also showed threshold-dependent sequence loss, including the loss of sequences repeatedly detected across replicate samples. Conclusions: Pipeline performance in highly diverse 16S rRNA amplicon datasets varies with dataset characteristics, evaluation criteria, and parameter settings. Pipeline selection should therefore reflect dataset characteristics and the analytical objective rather than rely on a single measure of performance, and minimum abundance thresholds should be evaluated in relation to sequencing depth rather than applied as fixed defaults. These findings emphasise the importance of reconsidering established analytical practices and transparently reporting bioinformatic settings as 16S rRNA amplicon sequencing is increasingly applied to highly diverse environmental microbial communities.

Passive sampling is the deployment of a collection material in the environment to continuously capture environmental DNA (eDNA) over time, offering the potential to integrate biodiversity signals while reducing the need for repeated active water collection. However, the mechanisms governing eDNA capture and retention on passive samplers remain poorly understood, limiting the interpretation of passive eDNA signals and their broader application. Here, we investigated the mechanistic performance of glass fibre passive samplers using controlled mesocosm experiments with three invasive freshwater bivalves: zebra mussels (Dreissena polymorpha), quagga mussels (Dreissena bugensis), and Asian clams (Corbicula fluminea). Specifically, we quantified eDNA accumulation dynamics, evaluated the contribution of different eDNA states, tested the persistence of captured eDNA, and compared passive sampler signals with conventional active sampling. Passive samplers rapidly accumulated target eDNA within hours of deployment, after which concentrations either plateaued or continued to increase depending on species. Sequential transfer of passive samplers between mesocosms containing different species showed that previously captured eDNA declined while new target eDNA accumulated to concentrations comparable to freshly deployed samplers, demonstrating continual turnover rather than permanent retention. Dissolved eDNA showed little evidence of accumulation beyond the concentration retained in the pore water within the membrane, suggesting that it is unlikely to be the dominant contributor to long-term passive sampler signals. Instead, the observed variability among replicate samplers, together with the physical properties of glass fibre membranes, suggests that membrane-bound and particulate eDNA are the primary contributors to passive eDNA capture. Collectively, these findings support a model in which glass fibre passive sampler signals reflect a dynamic equilibrium between ongoing eDNA capture and concurrent loss processes rather than cumulative accumulation over time. This mechanistic framework provides a foundation for interpreting passive eDNA data and informs the future development of passive sampling materials, deployment strategies, and biodiversity monitoring applications.

Environmental DNA (eDNA) metabarcoding enables sensitive, non-invasive assessment of fish communities, but highly multiplexed analyses using Oxford Nanopore Technologies (ONT) platforms require stringent control of sample-index misassignment and sequencing errors. We developed a MiFish experimental workflow combining unique dual indexes, a library protocol that omitted end preparation and used 5'-phosphorylated primers, BLAST-based demultiplexing, quality-dependent clustering, consensus generation, and haplotype partitioning by SNP/INDEL patterns. Omitting pooled end-prep reduced the mean index-chimera rate from 0.0674% to 0.000420%, a 160-fold reduction. We applied the workflow to an archive of seawater samples collected weekly off Monbetsu, Hokkaido, Japan, from 2012 to 2022. Relative read abundance (RRA) data were obtained for 274 samples, and eight taxa showed significant seasonality. For six of these taxa, the three-year mean RRA peaks coincided with reported spawning periods in Hokkaido. The workflow substantially reduced index misassignment and enabled cost-efficient, highly multiplexed analysis of a decadal fish eDNA time series.

Health monitoring is essential to ensure laboratory animal welfare and the reliability of experimental data in zebrafish research facilities. Conventional surveillance strategies based on resident or sentinel fish have limitations in terms of sensitivity and animal use, highlighting the need for alternative approaches consistent with the 3Rs principles. In this perspective, the present study evaluated the use of sludge collected from recirculating aquaculture systems as an environmental matrix for molecular health monitoring. Because sludge accumulates microorganisms and organic material from the entire system, it represents a promising sample for pathogen detection. Following an initial environmental surveillance phase to detect pathogens present in the system, the study aimed to optimise a molecular protocol for sludge analysis. To this end, four commercial DNA extraction kits were evaluated to assess their effectiveness in recovering bacterial DNA from sludge. Their performance was analysed by quantitative PCR in terms of extraction efficiency, repeatability, and limit of detection. The results highlighted the strengths and limitations of each DNA extraction protocol and confirmed the suitability of sludge as a non-lethal matrix for pathogen detection. These findings support the implementation of environmental monitoring as a practical and cost-effective alternative to sentinel-based surveillance, improving pathogen detection while reducing animal use in accordance with the 3Rs principles and Directive 2010/63/EU.

Detecting dugongs at low densities remains a major challenge for conservation monitoring, particularly in remote or turbid coastal habitats where conventional surveys are difficult. Environmental DNA (eDNA) offers a complementary approach, but robust species-specific tools are needed before it can be confidently applied to dugong surveillance. We developed and analytically validated a quantitative PCR TaqMan assay for detecting dugong (Dugong dugon) eDNA and evaluated its performance across controlled aquarium conditions in Australia and Abu Dhabi, and wild populations in northern Australia and the Arabian/Persian Gulf. The assay detected dugong eDNA in aquarium and field samples, with Sanger sequencing confirming amplification of dugong mitochondrial DNA from positive field detections. Our assay provides a robust and geographically transferable tool for targeted dugong eDNA detection, performing reliably across contrasting coastal environments and widely separated populations. This approach can strengthen local-scale surveillance of low-density or difficult-to-detect dugong populations.

1. Most species are neither discovered nor named, let alone included in analyses that require biological information such as trait measurements, images, ecological information and genome-scale data. Specimen-level DNA barcoding can help discover many of these species rapidly, but everything beyond discovery requires vouchers organized into putative species. Yet, existing barcoding workflows lack efficient techniques for voucher recovery, creating a post-barcoding bottleneck that limits the ability of converting barcoded specimens into biological knowledge. 2. Here we present a low-cost, open-source workflow consisting of two stages. The first safeguards barcoded specimens by separating them from DNA extracts and transferring them from microplates into ethanol-filled glass vials. The second converts the resulting voucher collection into a searchable physical resource by linking barcode-derived molecular Operational Taxonomic Unit (mOTU) assignments to vial positions and enabling specimens to be sorted into putative species either manually or automatically using a newly developed open-access robot (SORTER). 3. We evaluated the workflow using 2,024 insect specimens distributed across 21 96-well plates. For the first stage, DNA separation and specimen transfer required approximately 15 minutes per plate. For the second stage, MOTUmapper generated retrieval coordinates in a few seconds, after which the 2,024 vouchers belonging to the 452 putative species could be recovered manually in 5 days or with SORTER in 5 hours. Throughout both stages, specimen identities remained linked to barcode sequences, metadata and storage positions. 4. Vouchers are the Rosetta stones of biology because they connect different kinds of data to the same specimens. By safeguarding these vouchers and making them searchable, the workflow converts barcode projects from one-time molecular surveys into reusable resources for ecological and evolutionary research.

Hydropower development has emerged as a major driver of global environmental change, causing significant ecological damage to freshwater ecosystems and threatening aquatic biodiversity. However, the impact of cascade dams on fish communities in large river systems remains poorly understood, particularly in biodiversity hotspots such as China. Here, we used environmental DNA (eDNA) metabarcoding to assess the spatiotemporal patterns of fish diversity and community structure in the Ganjiang River, the largest tributary of the Yangtze River, in relation to cascade hydropower development. We collected water samples from 19 sites along the main stem and major tributaries in both wet and dry seasons, and analyzed the eDNA data using high-throughput sequencing and advanced bioinformatics pipelines. Our results revealed significant seasonal and spatial variations in fish diversity and community composition, with lower species richness and distinct assemblages in the impounded sections upstream of the dams compared to the downstream reaches. The fish communities were strongly influenced by environmental factors such as water temperature, pH, dissolved oxygen, and nutrients, which exhibited complex interactions with dam operation and land use practices. Spatial scale analyses further indicated that the impact of dams on fish diversity was most pronounced at the local scale, but also detectable at the regional scale. Our findings provide novel insights into the ecological consequences of cascade hydro-power development on fish communities and highlight the importance of considering spatial scales and environmental heterogeneity in assessing the impacts of human activities on freshwater biodiversity. Our study demonstrates the power of eDNA metabarcoding as a cost-effective and non-invasive tool for biodiversity assessment and conservation planning in large river systems facing multiple stressors and knowledge gaps.

Ranges of small benthic fauna are notoriously difficult to assess. In some of these cases, modern eDNA methods can shed light on species occurrence. Here we conduct an exploratory study on the fish eDNA recoverable from the gastrovascular cavities of the easy-to-sample pore water siphoning benthic invertebrate, Cassiopea, across six sites within the Florida Keys. Twenty-seven fish 12S identities were recovered from water samples, two from sediment samples, and seventeen from Cassiopea gut swabs. In total, thirty-two different species were identified from nineteen families, including one shark species (Ginglymostoma cirratum), and five species of cryptobenthic reef fishes (f: Gobiidae, Labrisomidae). Additionally, five species were identified from medusae samples that were not recovered in water or sediment samples. The species identities recovered may provide insight into the fish in direct proximity to Cassiopea assemblages, as well as indicate that Cassiopea may accrue disproportionate eDNA from cryptobenthic reef fish compared to surrounding environmental samples. The unorthodox sampling technique of using eDNA recovered from jellyfish stomachs yields another avenue for epibenthic community data acquisition.

Urbanization has been previously shown to affect evolution of plants and animals by restricting gene flow, increasing genetic drift, and causing divergent selection between rural and urban populations, however, the effect of urbanization on terrestrial gastropod evolution has only been minimally researched. Most research on terrestrial gastropods primarily focusing on invasive pest species, and little attention is given to native species like Ariolimax columbianus (Pacific banana slug). There is no current published research on the phylogeography of A. columbianus. In this study, surface epithelial cells from 66 Ariolimax columbianus individuals were collected from 11 locations around Corvallis, Oregon, and one site in northern California. These samples were amplified and sequenced at two loci: cytochrome oxidase 1 (CO1) and the internal transcribed spacer 2 (ITS-2). This data was used to generate maximum-likelihood phylogenetic trees to better understand the population genetics of these local native gastropods. The CO1 phylogenetic analysis showed two major haplotypes around the city of Corvallis - northwestern and southern, while the ITS-2 phylogenetic analysis showed three general haplotypes around the city - northern, western, and southern. The results in this study show the presence of genetic breaks between many of these populations, possibly due to a combination of factors including anthropogenic fragmentation, natural barriers, and the low active dispersal ability of A. columbianus. However, the extent that each of those factors has on disrupting gene flow and explaining patterns of genetic variation in this species requires further research. The present study acts as a foundation for the development of future research on banana slug phylogenetics and population genetic structure.

Batoids (rays, skates, and relatives) are among the most threatened cartilaginous fishes globally, yet the molecular diversity of Peruvian taxa remains poorly evaluated. Furthermore, regional barcode sequences for high-conservation-priority taxa, such as sawfish, manta, and devil rays, are often restricted to private records, hindering essential phylogeographic and population structure assessments. To bridge this gap, we developed BATOseq-PE, the first public genetic reference library for Peruvian marine batoids, covering 14 species sequenced in vitro from local specimens, complemented by a novel in silico assembled mitogenome for Notoraja martinezi. This open-access repository contributes 73 novel COI barcodes, providing robust taxonomic resolution across 15 molecular operational taxonomic units (MOTUs) and Barcode Index Numbers (BINs). Notably, we provide the first Peruvian public records for high-priority species such as Mobula munkiana, Mobula thurstoni, Myliobatis peruviana, Rostroraja velezi, and the first global public barcodes for N. martinezi and Zapteryx xyster. By coupling BATOseq-PE with global datasets, we conducted the first multi-taxa screening for delimitation and phylogeographic assessments on Peruvian batoids, revealing a cryptic lineage within R. velezi and genetic connectivity patterns suggesting female philopatry in Mobula mobular. Finally, by integrating nucleotide diversity into a Genetic Diversity Risk Indicator (GDRI) alongside ecological and fishery data, we identified an alarming trend where heavy fishing pressure overlaps with extinction risk and depleted genetic diversity, allowing us to suggest data-driven, tailored management strategies. BATOseq-PE represents a valuable resource to accelerate taxonomic and population assessments, ultimately fostering ray conservation through the integration of genetic frameworks into fishery management.

International wildlife trade represents one of the major drivers of biodiversity exploitation worldwide. However, the true taxonomic diversity embedded within commercial wildlife products often remains unknown because processing removes diagnostic morphological characteristics, preventing reliable species identification. Consequently, biodiversity assessments based solely on product labels may substantially underestimate the diversity of species involved in trade networks. To investigate hidden biodiversity within wildlife trade products, we applied DNA barcoding based on the mitochondrial cytochrome c oxidase subunit I (COI) gene to 77 products commercialized as fish swim bladders and seized at Guarulhos International Airport, Brazil. Molecular analyses successfully identified all samples and revealed the presence of four species: Plagioscion auratus (n = 38), Cynoscion acoupa (n = 7), Melanosuchus niger (n = 17), and Caiman crocodilus (n = 15). Fish species accounted for 71.4% of all samples, whereas crocodilians represented 28.6%, demonstrating that products marketed under a single commercial category may conceal substantial taxonomic diversity. Notably, the occurrence of two Amazonian crocodilian species within a trade chain traditionally associated with fish products reveals a previously undocumented component of the international wildlife trade. Our findings demonstrate that DNA barcoding is an effective tool for uncovering hidden biodiversity within processed wildlife products and provide evidence that wildlife trade networks may involve a broader spectrum of species than suggested by commercial labels. These results highlight the importance of molecular surveillance for biodiversity monitoring, wildlife trade regulation, and conservation planning.

Benthic algae integrate catchment-and reach-scale pressures in streams. We assessed land-use effects on benthic algal assemblage structure, chlorophyll-a, ash-free dry mass (AFDM), and autotrophic index across 24 headwater streams in the Sondu-Miriu River basin, Kenya, and explored wet-dry patterns in two monitoring streams. We sampled across four land-use categories: natural forest (NF), tea and tree plantations (TTP), smallholder agriculture (SHA), and smallholder tea (SHT). Synoptic sampling occurred in August 2025, and one NF stream and one SHA stream were sampled five times between November 2024 and May 2025. Shannon and Simpson diversity were highest in NF streams and lowest in SHT streams. Chlorophyll-a and AFDM were highest in SHA streams; chlorophyll-a was lowest in SHT, while AFDM and autotrophic index were lowest in NF. Relative abundance patterns showed that SHA streams were dominated by Lyngbya and Stigeoclonium, whereas SHT streams were dominated mainly by Stigeoclonium; NF streams supported persistent diatom-rich assemblages. Non-metric multidimensional scaling showed land-use separation of benthic algal assemblages, while redundancy analysis indicated weak, non-significant evidence of environmental associations in the synoptic dataset. Composition and biomass metrics tracked land-use gradients and wet-dry patterns, supporting periphyton-based bioassessment for monitoring disturbance in Afrotropical headwater streams.

The Community Conservation Index (CCI) developed by Chadd and Extence (2004) has been widely used to assess the conservation value of freshwater macroinvertebrate assemblages by integrating species rarity and community richness into a single metric. However, the original conservation scoring system was based largely on historical rarity designations and regional distribution data that no longer reflect current conservation assessments or recording coverage. This paper presents a comprehensive revision of the Conservation Scores underpinning the CCI. Species scores were updated using the latest IUCN-compliant Red List assessments for British freshwater invertebrates. For taxa classified as Least Concern, Data Deficient, or lacking a recent conservation assessment, scores were derived from national occurrence data across Great Britain for the period 1995-2024 using hectad occupancy. In addition, all non-native species were assigned a conservation score of zero to ensure they do not contribute positively to site conservation value. The revised dataset includes 1,680 taxa, compared with 1,117 in the original index, with 755 taxa receiving Conservation Scores for the first time. Thirty non-native taxa were assigned a score of zero. Among previously scored taxa, 165 received higher scores, 576 received lower scores, and 184 remained unchanged. The predominance of score reductions reflects improved knowledge of species distributions and the replacement of historical rarity classifications with contemporary assessments based on conservation status and national frequency of taxon occupancy. The revised CCI maintains the original strengths of combining rarity and richness while improving transparency, consistency and ecological relevance. By incorporating current Red List assessments, Great Britain-wide distribution data and revised treatment of non-native species, the updated index provides a more robust and nationally representative measure of freshwater macroinvertebrate conservation value. It offers an enhanced tool for biodiversity assessment, conservation prioritisation, site designation and environmental decision-making across freshwater ecosystems in Great Britain.

Freshwater biodiversity is declining at alarming rates globally, yet our capacity to monitor it depends heavily on the methods used for biodiversity assessment. Stream benthic macroinvertebrates are widely used as bioindicators, but traditional morphological identification is labor-intensive and limited by taxonomic expertise, while DNA metabarcoding offers a promising but incompletely evaluated alternative. Here, we compared biodiversity patterns detected by four macroinvertebrate sampling and identification methods (morphological identification of benthic samples, morphological identification of combined benthic and stream-edge samples, eDNA metabarcoding of bulk benthic samples, and eDNA metabarcoding of streamwater) across 62 sites in seven watersheds spanning the southern continental United States. We compared methods across three facets of biodiversity (gamma diversity, alpha diversity, and taxonomic composition) and three taxonomic levels (family, genus, and species). We found that which method detected more taxa shifted with taxonomic resolution: morphological methods detected more families and genera, while metabarcoding methods detected more species. Bulk benthic metabarcoding detected more families and genera than water eDNA metabarcoding, consistent with the "watered-down biodiversity" effect of DNA dilution in water samples, but the two metabarcoding methods did not differ significantly in species richness. All four methods detected largely distinct suites of taxa, with metabarcoding methods yielding more unique taxa than morphological methods. Notably, taxa uniquely detected by water eDNA were dominated by soft-bodied organisms (oligochaetes, leeches, earthworms) and water-column-associated taxa, while benthic eDNA uniquely detected hard-bodied EPT insects. We attribute this pattern to fundamental differences in eDNA shedding, transport, and deposition dynamics across organism types. These results demonstrate that no single method captures a complete picture of stream macroinvertebrate biodiversity, and that method choice should be guided by the taxonomic resolution and community components most relevant to the goals of the bioassessment program.

Full-length 16S rRNA gene sequencing using Oxford Nanopore Technologies has emerged as a promising approach to improve species-level resolution in microbiota studies. However, the accuracy of taxonomic assignment remains highly dependent on the bioinformatics s used to process Nanopore long-read data. Therefore, the only way to ensure a good level of certainty in obtained results is to use positive controls in the form of mock communities in the experimental designs. In this study, we compared the performance of Epi2Me 16S (using Minimap2 or Kraken2) workflows provided by Oxford Nanopore Technologies and an EMU workflow for full-length 16S rRNA gene analysis. Using a commercial mock community sequenced across multiple Nanopore runs, taxonomic assignment accuracy and reproducibility was evaluated. Epi2Me-Kraken2 exhibited 18 % of incorrect genus-level assignments and failed to identify 3 species present in the mock community. While Epi2Me-Minimap2 achieved an excellent genus-level classification, reporting 9 % of sequences assigned to a genus not in the mock community, species-level assignments were inconsistent for several community members such as Listeria. In contrast, EMU provided accurate and consistent species-level taxonomic profiles, with all species correctly identified while keeping the number of genus absent from the mock community at 1.2%. ImportanceThese results highlight that Epi2Me integrated workflows are not the best option for specie-level taxonomic assignation. More importantly, this paper underscores the importance of routine inclusion of positive controls for microbiota studies, in the form of mock communities, as a critical safeguard for accurate data interpretation. Without the use of a mock community, a paper published would be at risk of reporting wrong observations and inaccurate conclusions.